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I. INTRODUCTION 

In this work we employ a very recent addition to the range of moment methods available now to the problem 
of the s = 1/2 antiferromagnetic Heisenberg model on a sxiuare lattice. These other methods are based essentially 
on linked-cluster expansions and include the i-expansion |l|, the connected moment expansion (CMX) f^, and the 
coupled-cluster expansion 0, H]- The method we use is based on the Lanczos method, 



Hllpn) = Pn\lpn~l) +an|V^n) +/?«+! |'0«+l) 



(1) 



starting with a trial state |'0o)j whereby approximate yet analytic expressions for the Lanczos coefficients a„(-/V) and 
Pn{N) can be found at an arbitrary Lanczos step n in terms of low order cumulants, denoted by c„ n = 1, . . .. This 
expansion, also called the "plaquette expansion" |y , Q , is actually an expansion of the exact coefficients about the 
infinite size limit, that is in 1/A^. Furthermore in the combined limit of complete Lanczos convergence n ^ oo and 
the thermodynamic limit N ^ oo scaled Lanczos coefficients emerge in terms of z = n/N, taken to the lowest three 
orders 
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In Ref R] it was found that the ground-state energy density in the bulk limit to be given in general by 

eo = inf [a(z) — 2/3(z)] = inf e{z) . 



(2) 



(3) 



All these moment methods use the same input data, a sequence of connected moments up to some cut-off order, 
yet have rather different treatments of these so that it is of interest to compare them with each other. This work 
is the logical extension of Ref S which compared the t-expansion and the connected moments expansion with the 
coupled-cluster method, series expansions and spin-wave theory for the same model. In this work we have used the 
same set of cumulants as in the above work, up to the 15th order, and an additional set taken to the same order. 
However the interest of this work is not just confined to a discussion of the relative merits of the various methods, 
but also to demonstrate the contribution of moment methods to an understanding of the planar antiferromagnetic 
Heisenberg model. 



II. GROUND STATE AVERAGES FOR ANALYTIC LANCZOS EXPANSION 



Most previous applications of the plaquette expansion have been to the energy spectrum, notably the ground state 
energy ^ , |1^] , [^ and the mass gap ||l^ . Two earlier examples of the application of the expansion to other averages 
were the calculation of the staggered magnetisation of the isotropic 2D Heisenberg model |13] , Km . In these works 
the staggered magnetisation was extracted from the full nonlinear magnetic field dependence of the ground state 
energy for a system in a transverse external field. Here we take this approach to its logical conclusion and establish 
the formalism to describe an ground state average of an arbitrary operator V. Our approach is quite simple, we 
"piggy-back" our operator V onto our Hamiltonian, by tagging it with a variable A, and apply the earlier results and 
theorems to this new Hamiltonian. 

Consider the exact ground state energy of the new Hamiltonian H + XV , namely e(A) and with exact ground state 
wavefunction |^a}- 



ei\) = {^x\H + \V\-i'x) , 
then the Hellmann-Fcynman theorem gives the ground state average as 



(4) 
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So we need only evaluate the first order change in A in any quantity. Because we are assuming all quantities thus 
derived are analytic in some small neighbourhood of A = we require the overlap of the new ground state wavefunction 
with the original is nonzero, (^PoIvPa) 7^ 0. 

The new Hamiltonian H\ will generate a new sequence of Lanczos coefficients an{X) and /3^(A) from some suitable 
trial state |^o) and these in term will have analytic expansions 

an(A) = an + XS^an + 0{X^) 
PliX) ^Pl + XS^'PI + 0(A2) , 

where S^an and (5^/3^ define the first order shifts in a„ and /3^ in the presence of the operator V. In the interests of 
clarity a symbol without an argument of A is assumed to be the case where A = 0. Furthermore we consider models 
where the extensive scaling property when the Lanczos iteration number n and the system size N both tend to infinity 
with z = n/N fixed applies to the Lanczos coefficients, 

a«(A) "'^°° a{z) + Xd^a{z) + 0(A2) 
/3^(A)"-^-^°^/32(z) + A5VW + 0(A2). 

Then the exact ground state theorem ||7[] can be applied in the form 

e(A)=inf{«(A,z)-2[/32(A,z)]i/2} , (6) 

and one can show that 

{V) = viz)l^^, (7) 

where our new 2-function v{z) is 

v{z) = S a(z) — — , (8) 

and z is the 0-value of the minima in e{z) — a{z) — 2/3(z) if it exists. The new quantities in this, our central result, 
are the first order shifts in the Lanczos coefficients, and these will now be found in terms of moments. 

As is fundamental in the Lanczos process we need to find the moments of the new Hamiltonian with respect to our 
trial state, and the first order shifts of these T„, by 

{{H^D = (iJ") + Ar„ + 0(A2) , (9) 

and it is found that r„ is given by a sum of distributed generalised moments, 

n-l 

T„ = ^{H''-^-''VH^) , (10) 

fe=0 

when n > 1. The first member, Tq — 0. This result can be proven from the recurrence relation T„+i ~ (HTn) + {VH") 
or from the generating function 
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n=0 

From this point we go to the connected parts of these moments, with a first order shift denoted by Sn, 

{{H^r)c = {H'')c + A5„ + 0(A2) . (12) 

This connected part, or cumulant, has the standard generating function 
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and the standard recurrence relation with the moments 



n=l 



Tn+1 - 2^ ( ^ ](^ ')Sn+l-k 



(14) 



for n > 0. In effect what this reduces Sn to is a sum similar to Eq. |lO| except now it is over the connected parts of the 
distributed generalised moments, namely {H'^~'^VH^)c- It is these moments which are actually directly calculated. 
Finally, in the extensive many-body problem the connected moments have a size dependence via the coefhcients 



Sn = 5CnN 



(15) 



The expression for T„ (see Eq. ^^ or Sn in terms of the distributed generalised moments bears some similarity 
with the ones appearing in the t-expansion (see Equation 2.10 in Ref (§1), in that exactly the same averages occur, 
but that here the binomial coefficient is absent. This is because the Lanczos process is one where the Hamiltonian 
H generates new states in an enlarging state space, akin to a geometrical progression, and geometrical generating 
functions and convolutions arise. In contrast the i-expansion employs an exponentially mapped state e^^'^l-^p) and 
therefore exponential generating functions and convolutions appear. 

It should be noted that all the forgoing is exact, given that a system is solvable in this sense and all the moments, 
and Lanczos coefficients are analytically known. However for many non-trivial models one has knowledge only of the 
first say, 2r moments which give terms of z^ , . . . , z^~^ in a{z) and z, . . . ,z^ in 0^{z), or 2r -I- 1 moments which give 
terms in a{z) from z^ , . . . ,z^ and terms in 0^{z) from z, . . . , z^ . In the even case, with 2r moments, we define the 
truncation order as r and this is a natural ordering as there are equal numbers of terms in a and /3^ coefficients, and 
in the odd case 2r -I- 1 the truncation order is still r, but is not a natural order, but termed supplemented, in having 
one more term, namely a z^ in a{z). 

Given the shifted cumulant coefficients (5c„ , and therefore the shifted moments T„ there exists a number of algorithms 
for generating the shifted Lanczos coefficients, and one of the more robust and efficient ones is briefly described in 
the Appendix. The first few terms in the Lanczos coefficients are given by 
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Lastly we report some similar results to the above for the excited state gap, which was first found at the level of 
the first order plaquette expansion in Ref [|2| , and has been generalised in Ref. pj] . The connected moments for the 



excited (triplet) state c^N are related to the ground (singlet) state cumulants c^iV by 

One can define shifted Lanczos coefficients by 
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and there is an analogous exact result for the triplet gap 

9^rn{z)l=,. (19) 

where the gap function m(z) is defined as 

m{z) = 5^a{z) - '-^^ . (20) 

The shifted cumulants 5'^Cn can be found from the equivalent i-expansion function R{t) B via 



III. APPLICATION TO THE 2D XXZ HEISENBERG MODEL 

We apply the plaquette expansion formalism to the s = 1/2 anisotropic antiferromagnetic Heisenberg Model on a 
square lattice, with the usual Hamiltonian 



11 — y o^ o j n- X\ o^ o j ~r Oj^ o 



(22) 



<«j> 



where x — Q denotes the Ising Model and x = 1 the isotropic antiferromagnetic Heisenberg Model. For the trial ground 
state we take the "natural" and commonly employed choice of the eigenstate to a; = problem, the classical Neel 
state, and for the excited triplet trial state a state with one spin of the Neel state flipped. We have used the cumulants 
/„ and Rn in Ref . g , and generated a new set of staggered magnetisation cumulants to an order equivalent to that 
used in the above reference. From these we have generated the Lanczos coefficients and the staggered magnetisation 
Lanczos coefficients up to 7th order, and the triplet gap Lanczos coefficients up to 5th order. All new cumulants and 
the Lanczos coefficients are displayed in the Appendix. 

Using these cumulants the ground state energy density eo, the staggered magnetisation A4 and triplet gap Q have 
been computed for a set of anisotropy parameters and for all truncation orders up to the maximum, by evaluation at 
the minima of e{z). This is the simplest way of analysing the plaquette expansion and is free from any biasing and 
assumptions, so it is the preferred option. Given more understanding of the convergence properties of the plaquette 
expansion a more sophisticated analysis or extrapolation strategy may be employed, and we refer the reader to Ref 
|ll| for a discussion of this issue. The results for the ground state energy density, the staggered magnetisation and the 
excited state gap are displayed in Table |, Table ||, and Table [II respectively. Depending on the order and anisotropy 



parameter, the minima in the ground state energy function e[z), could be complex, usually close to the positive real 
axis, however. In these cases we took the real part of the function computed at the complex minima, and these are 
marked with a asterix(*). To treat the case of no real minima properly some form of extrapolation would be needed 
but we do not pursue this issue here for several reasons. Firstly, a real minima exists for all the data at the highest 
order (r = 7), which is the most important, at least for the ground state energy and magnetisation. Secondly there is 
some continuity at a given order, between the data where a real minima exists and that data where it is complex. As 
the anisotropy parameter varies through these regions the minima shifts off slightly from the real axis but nonetheless 
remains close. There is also a second number appended to each data entry of the tables, and that is the difference 
between that value computed using the natural ordering and the supplemented ordering. This is the most unbiased 
estimate we have for indicating the systematic error in our expansion although there are no rigorous results to say 
that this difference somehow bounds the error - it is at best semi-quantitative. 

The plaquette expansion results display two trends - one where as the truncation order increases we find substantial 
and systematic improvement in the averages, at least over the region where a real minima exists, and another as 
the anisotropy parameter increases from zero, where we find rapidly increasing errors, although the extent of this 
is dependent on the particular quantity. The groimd state energy density is the most accurate, then the staggered 
magnetisation, and the triplet gap is the worst. 



IV. COMPARISON OF THE METHODS 



Taking the highest order plaquette expansion results, we have tabulated these against the other methods in Table |IV|, 
Table and Table 



For the ground state energy one can see clear trends in the comparison : 

• for a: < 0.5 the i-expansion and CMX are slightly better than the plaquette expansion although the actual 
differences between them and with the series results are very small, 

• for X — 0.8 the plaquette expansion is better than the t-expansion and the CMX, while for x — 0.9 the i-expansion 
is a slight improvement over the plaquette expansion, 

• and for x > 0.95 the plaquette expansion is clearly better than the f-expansion, the CMX and 3rd order spin- wave 
theory, and especially so at a; = 1. 

For the staggered magnetisation one observes the following trends - 

• for X = 0.5 the plaquette expansion is better than the i-expansion and the CMX, 

• for X = 0.8 this ordering has reversed with the i-expansion superior to the plaquette expansion and the CMX, 
while for x — 0.9 the order gets reversed once more, and the plaquette expansion is a slight improvement over 
the i-expansion, 

• and for a; — > 1.0 the plaquette expansion stabilises on a higher value than that given by the t-expansion, 3rd 
order spin-wave theory, and the series expansion, although remaining less than the CMX result. 

And for the triplet gap one finds that all moment methods give a rather high nonzero value at the isotropic point, 
whereas its vanishing for the spin-wave theory and series expansion is built-in. We would like to emphasise that for 
X > 0.9 there was no real minima at the highest order that applied in the plaquette expansion, and consequently the 
estimated errors grossly underestimate the true error. The trends here are - 

• for a; = 0.2 — > 0.5 the i-expansion is better than the CMX and the plaquette expansion, 

• for X = 0.8 —>■ 0.9 the f-expansion is still superior to the plaquette expansion although it, in turn is better than 
the CMX, 

• and for x > 0.95 the plaquette expansion is slightly better than the other moment methods but all have errors 
which swamp the actual values. 

Finally we compare these results with the coupled-cluster results at the isotropic point. The results given in 
Table IV and Table |^ are the extrapolated results based on the first few orders of the approximation, and differ 



significantly from the highest order results. Both recent extrapolated ground state energies are better than either the 
i-expansion or the CMX and comparable in accuracy to the plaquette expansion energy. The extrapolated staggered 
magnetisations are also very close to that predicted by the plaquette expansion, probably higher than the true value. 

The relatively poor performance of the i-expansion arises from a number of reasons. Firstly the extrapolation t —^ oo 
is done without any knowledge of the global analytic properties of E(t), and it is widely known in extrapolation 
work that this is perhaps the largest single source of poor convergence and inaccuracies. In contrast there is no 
extrapolation problem in the plaquette expansion, so long as a real minima exists. If this is not the case then some 
form of extrapolation of the kind employed in Ref pT[ may be necessary here. 

Another source of mediocre performance, and this is shared with the plaquette expansion, is the rather poor quality 
of trial state. The classical Neel state is usually chosen because it is simple and relatively straightforward to generate 
moments with, but a trial state which can generate the correct singularity structure at the isotropic point may give 
better results even though fewer moments may be found from it. 

In summary we found good quantitative predictions from the plaquette expansion for this model but, as expected, 
the worst results occured near the isotropy point where a first order transition to a gapless state is expected. The 
plaquette expansion was clearly better at predicting the ground state energy than the other moment methods, had a 
slighter higher staggered magnetisation than the others, and all of the methods gave a nonzero gap at the isotropic 
point. We would like to emphasise that many aspects of the moment methods are still poorly understood, and that 
considerable scope exists for improving their accuracy. 

We would like to acknowledge that this work was supported by the Australian Research Council. 



V. APPENDIX 

Using efficient graph enumeration and calculation of the embedding constants Ref. |q| , [16| we have the shift in the 
staggered magnetisation cumulants 
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and using the excited state coefficients from Ref. pi we have the shift in the excited state cumulants 
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The expansion for the Lanczos coefficients are given exactly ||l^ , [Q up to their truncation order by 
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The expansion for the change in the Lanczos coefficients corresponding to the staggered magnetisation is given by 
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The expansion for the shift in the Lanczos coefficients describing the triplet energy gap is given by 
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TABLE I. The ground state energy eo and estimated error for the 2D XXZ Model calculated using the plaquette expansion 
as a function of the anisotropy parameter x and the order of truncation r. Those energies where one has the case of a real 
minima are unmarked, whereas those energies which are given by the real part evaluated at a complex minima are marked by 
a asterix. 
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-0.60693(2) 


0.9 


-0.622(8) 


-0.633(2) 


-0.6363(3) * 


-0.63667(7) * 


-0.6365(2) * 


-0.6360(1) 


0.95 


-0.64(1) 


-0.648(4) 


-0.6526(2) * 


-0.65280(3) * 


-0.653(1) * 


-0.6521(4) 


0.98 


-0.64(1) 


-0.658(4) 


-0.6627(1) * 


-0.66276(2) * 


-0.663(1) * 


-0.6622(5) 


0.99 


-0.65(1) 


-0.661(5) 


-0.66604(9) * 


-0.66613(4) * 


-0.666(1) * 


-0.6657(6) 


1.0 


-0.65(1) 


-0.664(5) 


-0.66945(7) * 


-0.66952(7) * 


-0.670(1) * 


-0.6691(6) 
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TABLE II. The staggered magnetisation M and estimated error for the 2D XXZ Model as calculated using the plaquette 
expansion as a function of the anisotropy parameter x and the order of truncation r. Again those energies where one has the 
case of a real minima are unmarked, whereas those energies which are given by the real part evaluated at a complex minima 
are marked by a asterix. 



M 


r = 2 


3 


4 


5 


6 


7 


a; = 0.2 


0.49563(7) 


0.495543(8) 


0.495532(2) 


0.4955293(8) 


0.4955281(4) 


0.4955275(2) 


0.5 


0.475(2) 


0.4717(4) 


0.47112(8) 


0.47100(2) 


0.470960(9) 


0.470946(4) 


0.8 


0.44(1) 


0.424(6) 


0.416(1) 


0.4143(4) 


0.4152(6) 


0.4162(4) 


0.9 


0.43(2) 


0.40(2) 


0.373(5) * 


0.3794(4) * 


0.377(3) * 


0.383(3) 


0.95 


0.43(2) 


0.39(3) 


0.367(4) * 


0.371(2) * 


0.37(1) * 


0.366(6) 


0.98 


0.42(3) 


0.38(2) 


0.364(3) * 


0.367(3) * 


0.359(1) * 


0.358(5) 


0.99 


0.42(3) 


0.37(2) 


0.363(3) * 


0.365(3) * 


0.35676(7) * 


0.355(5) 


1.0 


0.42(3) 


0.37(2) 


0.362(3) * 


0.364(4) * 


0.355(4) * 


0.353(5) 
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TABLE III. The excited state gap Q and estimated error for the 2D XXZ Model calculated using the plaquette expansion 
as a function of the anisotropy parameter x and the order of truncation r. Again those energies where one has the case of a 
real minima are unmarked, whereas those energies which are given by the real part evaluated at a complex minima are marked 
by a asterix. 



r = 2 



a; = 0.2 

0.5 

0.8 

0.9 
0.95 
0.98 
0.99 

1.0 



1.943(4) 
1.67(4) 
1.2(1) 
1.1(2) 
1.0(2) 
1.0(2) 
1.0(3) 
1.0(3) 



1.938(1) 
1.617(7) 
1.03(5) 
0.8(2) 
0.6(2) 
0.5(2) 
0.5(1) 
0.5(1) 



1.936(5) 
1.607(2) 
0.962(9) 
0.56(5) * 
0.47(3) * 
0.42(3) * 
0.40(2) * 
0.39(2) * 



1.936 
1.603 
0.953 
0.61 * 
0.51 * 
0.45 * 
0.42 * 
0.40 * 
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TABLE IV. Comparison of the most accurate plaquette expansion values for the ground state energy eo with those of the 
i-expansion, the connected moments expansion, spin-wave theory and series expansions, taken for various anisotropy parameters 

X. 



£0 


plaquette 


i-expansion^ 


i-expansion'' 


CMX" 


3rd Order'* 


series" 


coupled' 




expansion 


D Fade 


Laplace 




spin-wave 


expansion 


cluster 


a; = 0.2 


-0.5066651(1) 


-0.5066653 


-0.5066653 


-0.50666529 


-0.50657179 


-0.5066653 




0.5 


-0.541631(1) 


-0.541636(3) 


-0.5416359 


-0.54163641 


-0.5413803 


-0.5416371 




0.8 


-0.60693(2) 


-0.6068(2) 


-0.6067604 


-0.60677223 


-0.607376 


-0.606902(2) 




0.9 


-0.6360(2) 


-0.6357(2) 


-0.6353801 


-0.63537633 


-0.636654 


-0.635844(4) 




0.95 


-0.6521(2) 


-0.6518(3) 


-0.6510589 


-0.65101764 


-0.652718 


-0.65189(1) 




0.98 


-0.6622(2) 


-0.6618(4) 


-0.6609227 


-0.66083842 


-0.66287 


-0.66211(2) 




0.99 


-0.6656(3) 


-0.6651(6) 


-0.6642888 


-0.66418527 


-0.66637 


-0.66563(6) 




1.0 


-0.6691(3) 


-0.668(1) 


-0.6676946 


-0.66756890 


-0.6699993 


-0.6693(1) 


-0.6692 
-0.6691(3) 



''Reference 
''Reference 
'^Reference | S| 
"^Reference ic| 
•^Reference |2C 
Reference [21 



TABLE V. Comparison of the most accurate plaquette expansion values for the staggered magnetisation M with those 
of the i-expansion, the connected moments expansion, spin-wave theory and series expansions, taken for various anisotropy 
parameters x. 



M 


plaquette 


i-expansion^ 


CMX"" 


3rd Order'^ 


series'* 


coupled'^ 




expansion 


D Fade 




spin-wave 


expansion 


cluster 


a; = 0.2 


0.4955275(2) 


0.495527(1) 


0.49552690 


0.49573699 


0.4955265 




0.5 


0.470946(4) 


0.4710(4) 


0.47097192 


0.47172243 


0.4709287 




0.8 


0.4162(4) 


0.4173(6) 


0.41727472 


0.416390 


0.416896(5) 




0.9 


0.383(3) 


0.386(4) 


0.39065099 


0.383864 


0.38553(2) 




0.95 


0.365(3) 


0.36(1) 


0.37519399 


0.3607157 


0.36266(6) 




0.98 


0.357(3) 


0.35(2) 


0.36517718 


0.340646 


0.3422(2) 




0.99 


0.355(4) 


0.34(2) 


0.36170935 


0.33068 


0.3319(4) 




1.0 


0.353(5) 


0.33(3) 


0.35817561 


0.3069 


0.307(1) 


0.35 
0.340(5) 



^Reference 
'^Reference 
'^Reference 
"* Reference 
"Reference 
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TABLE VI. Comparison of the most accurate plaquette expansion values for the excited state gap Q with those of the 
i-expansion, the connected moments expansion, spin-wave theory and series expansions, taken for various anisotropy parameters 

X. 



Q 


plaquette 


t-expansion^ 


CMX'' 


3rd Order'^ 


series'^ 




expansion 


D Fade 




spin-wave 


expansion 


a; = 0.2 


1.9357(6) 


1.9338(2) 


1.93383753 


1.942248 


1.933815 


0.5 


1.603(3) 


1.594(10) 


1.59431537 


1.629782 


1.59736(4) 


0.8 


0.953(9) 


0.96(2) 


0.99272985 


0.98798172 


0.970(3) 


0.9 


0.61(4)* 


0.65(5) 


0.74000985 


0.65721412 


0.66(1) 


0.95 


0.50(3)* 


0.52(10) 


0.60488156 


0.44092456 


0.45(3) 


0.98 


0.44(3)* 


0.45(15) 


0.52114736 


0.26588297 


0.26(6) 


0.99 


0.42(5)* 


0.4(2) 


0.49280432 


0.18381204 


0.17(8) 


1.0 


0.40(5)* 


0.3(3) 


0.46424870 


0.0 


0.0(1) 


''Reference 


1 










'^Reference 


i 










'^Reference | 


S 










"^Reference 


20| 
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